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ABSTRACT 

We present template radial velocity curves of afe-type RR Lyrae stars con- 
structed from high-precision measurements of Ha, H/3, and H7 lines. Amplitude 
correlations between the Balmer line velocity curves, Johnson \^-band, and SDSS 
g- and r-band light curves are also derived. Compared to previous methods, these 
templates and derived correlations reduce the uncertainty in measured systemic 
(center-of-mass) velocities of RR Lyrae stars by up to 15 km s _1 , and will be of 
particular interest to wide-area spectroscopic surveys such as the Sloan Digital 
Sky Survey (SDSS) and LAMOST Experiment for Galactic Understanding and 
Exploration (LEGUE). 

Subject headings: stars: horizontal-branch - stars: Population II - stars: vari- 
ables: RR Lyrae - techniques: radial velocities 



1. Introduction 



RR Lyrae stars are old (> 9 Gyr), low-mass (~ O.7M ), pulsating stars (V-band am- 
plitudes of Ay ~ 1 mag and p eriods of ~ 0.6 days) that reside in the instability strip of the 
horizontal branch ( )Smithlll995l ). They are valuable tracers of old stellar populations and have 
recently been used to map the Galactic halo structure (i.e., its stellar number density profile) 
and substructure (e.g., tidal streams and halo overdensities; 



20081 . iMiceli et al.ll2008l . IWatkins et al.ll2009l . ISesar et al.ll2010af ). 



Vivas et al. 2001, Keller et al. 



While the spatial distrib ution of RR Lyrae s tars is a powerful indicator of candidate halo 
substructures (e.g., Fig. 11 of lSesar et al.ll2010a[ ). it does not always reveal the full nature of 
halo substructures. For example, the Pisces Overdensity, a spatial group of about a dozen 



RR L yrae stars at ~ 80 kpc from the Sun (clump "J" in ISesar et al. 



20091 ). was initially suspected to be a part of a tidal stream ( ISesar et al.ll2010al ). However 



2007 



Watkins et al. 
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spectroscopic observations by iKoUmeier et al.l (120091 ) and ISesar et al.l ( 12 01 Obi ) have revealed 
that the overdensity actually consists of two velocity groups; one moving towards and the 
other one moving away from the Galactic center. This result suggests that the Pisces Over- 
density may be a pile-up of tidal debris near the orbital turn-ar ound point of a now -disrupted 
progenitor (which was, most likely, a dwarf spheroidal galaxy; ISesar et al.ll2010bl ). 



The IKoUmeier et al.l (120091 ) and ISesar et al.l ( 12 01 Obi ) findings underline the importance 
of obtaining kinematic data on candidate halo substructures; the nature and interpretation of 
a halo substructure can quickly change once kinematic data become available. Yet, relative 
to photometric data, which is abundant in this era of wide-area and multi-epoch surveys, 
spectroscopic observations are more scarce and may be considered as more "expensive" in 
terms of telescope time. Obtaining the kinematic information for variable stars, such as RR 
Lyrae stars, is even more complicated as one needs to subtract the velocity due to pulsations 
to get the center-of-mass (hereafter, systemic) velocity. 

There are two approaches to this problem of determining the systemic velocity of RR 
Lyrae stars. In the first approach, one tries to observe RR Lyrae stars at a particular point 
during a pulsation period when the observed radial velocity is presumed to be equal to the 
star's systemic velocity (usually at phase ~ 0.5 ± 0.1, see Section [6]). Unfortunately, some- 
times it is very difficult to schedule observations to target this point (e.g., when executing 
observations in queue or service mode). In addition, RR Lyrae stars are dimmer near phase 
0.5 (~ 0.3 — 0.5 mag compared to maximum light at phase of 0), and thus require longer 
exposures to reach the same signal-to-noise ratio. The danger of longer exposures is blurring 
of spectral lines due to pulsations, which in turn decreases the precision of radial velocity 
measurements. 

In the second approach, one assumes a certain model (a template) that describes the 
radial velocity due to pulsations as a function of pulsation phase. The amplitude, A rv , of 
the radial velocity curve template (hereafter, template) is then scaled using a correlation 
between the amplitudes of velocity curves and light curves. The velocity at any phase 
can then be easily read-out once the scaled template is matched to one or several radial 
velocity observations taken at known phases. An example of such a tem plate and an equation 
relating velocity and light curve amplitudes is provided by iLiul ( Il99ll ) (see his Equation 1). 
The advantage of this approach is that the systemic velocity can be estimated using radial 
velocities obtained in phases different from phase of 0.5. 

Unfortunately, Liu's template and his velocity vs. light curve amplitude relation are 
not really suitable for estimating systemic velocities of distant (i.e., faint) RR Lyrae stars 
because they were derived from measurements of metallic lines and should only be used 
with measurements of such lines. This is a problem since measuring radial velocities using 
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metallic lines requires high-resolution spectroscopy (R > 10000), and obtaining quality, 
high-resolution spectra of faint (V > 18) stars with exposures that are a small fraction of 
the pulsation period (< 10% to avoid spectral blurring) is virtually impossible with today's 
facilities. 

Instead of using weak metallic lines, a more practical choice would be to use strong 
Balmer lines (Ha to H<5), which can be observed with more widely available low- resolution 
(R ~ 1500) spectrographs. However, Liu's template and his velocity vs. light curve amplitude 
relation cannot be used with measurements of Balmer lines because the Balmer line vel ocity 
curve s have a different shape and amplitude than t he velocity curve of met a llic lines ( Qke 



1962 ]) . A metho d adopted by ma n y studies (e.g.. lHawlev &; Barnes! Il985l ; iLaydenl Il994 



Vivas et al. 2005; Prior et al 



2009 



curve of RR Lyrae star X Ari measured from the H7 line by lOkd (119661 ). This is not an ideal 



Sesar et al.l l2010bl . 120121) is to us e the radial velocity 



solution b ecause the velocity curve amplitude of X Ari (A 

rv ~ 85 km s , see Figure 1 of 



Okdll966l ) will likely not be similar to the velocity curve amplitude of some other star with a 
different V^-band light curve amplitude (Ay). The consequence of using an inadequate model 
for radial velocities is a more uncertain estimate of the systemic velocity. 

The ability to precisely measure systemic velocities of RR Lyrae stars is very important 
as it may help with various Galactic studies, such as identifying halo substructures using 
kinematics (e.g., the Pisces Overdensity groups), discerning the nature of progenitors of halo 
substructures (e.g., globular cluster vs. a dwarf sp heroidal galaxy) , and constraining orbits 



of ha lo substructures (e.g., the Virgo Overdensity; ICasetti-Dinescu et al.ll2009l ; ICarlin et al. 



20121 ). In addition, obtaining precise systemic velocities of RR Lyrae stars using one or two 
measurements would be beneficial to spectroscopic sur veys such as the LA MOST Experiment 
for Galactic Understanding and Exploration (LEGUE; lDeng et al.ll2012l ). which may observe 
several thousands of RR Lyrae stars within 30 kpc of the Sun, and which plan to study the 
kinematics and structure of the Galactic halo. 

In this paper, we use high-precision measurements of Ha, H/3, and H7 lines of six field 
afr-type RR Lyrae stars to first construct template radial velocity curves (Section [3]). We 
then derive correlations between amplitudes of Balmer line velocity curves and light curves 
observed in Johnson \^-band (Section Hj) or SDSS g- and r-bands (Section ED- The systemic 
velocity of RR Lyrae stars and its uncertainty are defined in Section [BJ and the results, along 
with implications for surveys such as LINEAR, PTF, LEGUE, and LSST, are discussed in 
Section [7J 
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Data set 



The sample of stars used in the construction of velocity templates consists of 6 field 
RR Lyrae type ab stars (see Table [I]) that were observed between 2006 and 2009 by George 
W. Preston with the echelle spectrograph of the du Pont 2.5 m telescope at the Las Campanas 
Observatory. The stars have l/-band light curve amplitudes (Ay) ranging from 0.64 mag 
(Z Mic) to 1.13 mag (RV Oct ), as measured by th e All Sky Automated Survey (AS AS; 
Pojmahskilll997l ) . According to ISzczygiel fc Fabry ckv (20 071) , three of the st ars in this sample 



exhibit the Blazhko effect in their light curves ( jBlazhko 



19071 ISmithlll995h . 



The observations_j_data reduction, and measurement of radial velocities are described in 



detail in iFor et al.l ( 120111 ). The data, kindly provided to us by George W. Preston (private 
communication), consist of radial velocity measurements obtained using metallic and Balmer 
lines (Ha, H/3, and H/y). The formal random error is ~ 0.5 km s _1 for radial velocities 
obtained using metallic lines, and ~ 3 km s _1 for radial velocities obtained using Balmer 
lines. In addition, systematic errors of ~ 1 to 2 km s _1 may also be present. The phase of 
each observati o n was calculated using periods and epochs of maximum light listed in Table 
5 of IFor et al.l (120111 ). with the time of maximum light corresponding to phase of 0. The 
number of observations per star ranges from about 100 (DT Hya) to about 240 (XZ Aps). 



3. Construction of velocity templates 

We start by normalizing radial velocities to the range [0, 1]. The normalization is done 
on a per star basis using 

TV f "V m i n , » 

T"V ma x Wmin 

where rv min and rv max are, respectively, the lowest and highest radial velocity determined 
after discarding the top and bottom 1% of radial velocities of that star. The normalized veloc- 
ities from all stars are then combined and a cubic B-spline is interpolated through the data. 
The normalized, phased, and combined radial velocities are shown in Figure [1] (top plots in 
each panel). The velocities and the interpolated B-spline (the template) are offset so the 
velocity at phase 0.5 (the phase at which the radial velocity is presumed to be equal to the sys- 
temic velocity) is equal to 0. The templates are provided in the supplementary data available 



in the electronic edition of the journal and at http : //www. astro . caltech. edu/~bsesar/RVtemplates . taj 



Note that the templates of Balmer line velocity curves do not cover the full range of 
phases, but end at phase of 0.95. Due to a rapid change in velocities of Balmer lines between 
phases of 0.95 and 1.0, we were not able to properly model the entire velocity curve with a 
single cubic spline. We do not consider this to be a major problem from an observational 
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point of view because the velocities measured from Balmer lines between phases of 0.9 and 
1.0 are in general more uncertain, mainly due to increased spectral blurring (i.e., due to a 
rapid change in velocity as a function of time), and are not very useful. 

The interpolated cubic B-spline defines the average shape of the radial velocity curve. To 
estimate the uncertainty in this average shape, we calculate the root-mean-square (hereafter, 
rms) scatter around the template and plot its dependence on phase in bottom plots in each 
panel of Figure [TJ The shape of the Ha template is most uncertain, with an rms scatter of 
~ 10% of A rv at phase of about 0.55 (~ 10 km s -1 for a RR Lyrae star with A v = 1.0 mag). 
The shape of the metallic lines template curve is the least uncertain (~ 3% of A rv ), followed 
by H7 and H/3 template curves (~ 3% to 4% of A rv or ~ 4 km s -1 for a RR Lyrae star with 
A v = 1.0 mag). 



4. A rv vs. Ay correlations 

Having defined average velocity curves (templates), we now fit templates to observed 
radial velocities of each star to determine the amplitudes of metallic and Balmer lines velocity 
curves (A rv ). The best-fit amplitudes are listed in Tabled! and are plotted versus l/-band 
amplitudes in Figure HJ 



As iLiul ( 119911 ) (hereafter, L91) first noted, there is a tight correlation between the 
amplitudes of metallic lines velocity curves and V^-band light curves. We also see a similar 
correlation, but for velocity curves measured from Balmer lines. Following L91, we do an 
unweighted, linear least-squares fit to velocity amplitudes as a function of the V^-band light 
curve amplitudes and find: 



A™ et = 25.6(±2.5)A y + 35.0(±2.3), a fit = 2.4 km s' 1 (2) 

A* 01 = 35.6(±2.5)A y + 78.2(±2.4), a fit = 3.4 km s' 1 (3) 

A^f = 42.1(±2.5)A V + 51.1(±2.4), a fit = 3.0 km s' 1 (4) 

A^ = 46.1(±2.5)A y + 38.5(±2.4), a fit = 2.8 km s~\ (5) 

where crgt is the rms scatter of the best fit. 

In the top left panel of Figure [2J we compare our best-fit for metallic lines (Equation [2]) 

to 

L91 _ 40.5A y + 42.7 
Arv - ^ > W 

which is the best-fit for metallic lines obtained by L91 (his Equation 1) scaled by the so- 
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called "projection factor" p = 1.37. Equation 1 of L91 needs to be scaled because L91 uses 
pulsation velocities, which are related to observed radial velocities as t> f>s = Vpuis/p (see Liu's 
Equation 2). The projection factor was estimated as the average ratio of velocities calculated 
using Liu's Equation 1 and metallic A rv values listed in T able U The es timated projection 



factor p = 1.37 is in the range of values listed by L91 and iKovacd (120031 ) (1.31 to 1.37). In 
conclusion, the scaled L91 fit and our fit are remarkably similar (within ~ 2.5 km s _1 of 
rms scatter), even though Liu's sample has more stars and spans a much greater range of 
y-band amplitudes than our sample (Liu's sample extends up to A v ~ 1.3 mag). 

Figure [2] shows that for a fixed l/-band amplitude, each of the Balmer lines studied 
here has a different velocity curve amplitude. This suggests that, in theory, the most precise 
systemic velocity may be obtained by averaging systemic velocities which are estimated from 
radial velocities measured from individual Balmer lines. However, in some cases (e.g., for 
faint stars), it may make more sense to measure radial velocities by cross-correlating more 
than one Balmer line with a reference spectrum. Due to the fact that Balmer line velocity 
curves have different shapes and amplitudes, we expect that such measurements may be 
more uncertain. 

To estimate the level of uncertainty introduced by using more than one Balmer line 
in a cross-correlation, in Figure [3] we show the standard deviation of radial velocities of 
various Balmer lines as a function of phase for RR Lyrae stars with V^-band amplitudes 
Ay = 0.6, 1.0, 1.4 mag. In general, the uncertainty in velocity introduced by using more 
than one Balmer line in a cross-correlation is lowest for phases earlier than 0.6 (< 4 km s -1 
for Ay = 1 mag), and it decreases with increasing V^-band amplitude. 



5. A K and A r vs. A v correlations 



Equations [2] to [5] describe the amplitude correlations between radial velocity curves and 
the l/-band light curve. To enable use of these eq uations with RR L yrae stars that have been 
observed in the Sloan Digital Sky Survey (SDSS; lYork et al.ll2000l ) g- and r-band filters, we 
also derive the correlations between the light curve amplitudes of RR Lyrae stars observed 
in SDSS g and r bands and the Johnson V band. 

In deriving the correlations, we use the best-fit templates and light curve parameters o f 
379 a6-type RR Lyrae stars observed in SDSS Stripe 82 (see Table 2 in lSesar et al.ll2010al ). 
We synthesize each RR Lyrae star's Johnson ^/-band light curve from the best-fit SDSS g- 
and r-band light curves (gLc(4>) an d rLC (</>)) using a photometric transformation given by 
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Equation 10 from llvezic et al.l (120071 ): 



Vlc(^) = 9lc{^) + O.O688^r LC (0) 3 - 0. 2056 gr LC (<\>) 2 - 0.3838^(0) - 0.053, (7) 

where gvLcifi) = 9lc((P) ~ r LC*(0), an d [0, 1) is the phase of pulsation. The V^-band 
amplitude is then simply measured from the synthesized light curve. 

Having measured V^-band amplitudes (Ay), we can now study how they correlate with 
SDSS g- and r-band amplitudes. We find that, within ~ 0.02 mag of uncertainty (rms 
scatter), the \^-band amplitudes can be calculated as 0.9A g and 1.2lA r , where A g and A r 
are SDSS g- and r-band light curve amplitudes of RR Lyrae stars. The 1.21 factor is also 
obtained for the Johnson R band, by comparing Johnson V- and -R-band light curves of RR 
Lyrae stars V3, V4, V5, and V10 observed in NGC 5053 bv lFerro et all J201oh . 



6. Systemic velocity and its uncertainty 

In Section [31 the templates were defined to have a zero point at phase of 0.5, since we 
assumed the radial velocity at phase of 0.5 to be equal to the systemic velocity. Using this 
definition, the observed radial velocity v D b s can be calculated as 

v bs(®obs) = A rv T(® obs ) + v lt (8) 

where A rv is the amplitude of the radial velocity curve (calculated using Equations [2] to [5]), 
T( < l ) obs) is the template radial velocity curve, $ Q b s is the phase of observation, and v 7 is the 
systemic velocity. 

We use the following equation to model the uncertainty in the systemic velocity (c y ) 

a Vi = (A rv + a flt ) 2 [a 2 emplate (^ obs ) + (0.1k)% (9) 

where k is the slope of a template between phases 0.4 and 0.6 (1.16, 1.47, 1.54, and 1.42 
for the metallic, Ha, H/3, and H7 line templates, respectively). In the above equation, the 
a tem P iate(^obs) term describes the rms scatter in the template at the phase of observation 
(bottom plots in Figure [1]), and the (0.1k) 2 term describes the error introduced by the 
uncertainty in phase at which the radial velocity is equal to the systemic (center-of-mass) 
velocity. These two terms are given in units of velocity amplitude, so they scale with A rv . 
To take into account the uncertainty in Equations [2] to O we add the rms scatter of the A rv 
vs. A v fit to our estimate of the velocity amplitude. 

The (0.1k) 2 term describes our lack of knowledge of the exact phase at which the radial 
velocity is equal to the systemic (center-of-mass) velocity. In Section [31 we assumed this 
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phase to be at 0.5. However, lOkd (119621 ) used radial velocity measurements of RR Lyrae 
star SU Dra and found the systemic velocity of this star to be consistent with the observed 
radial velocity at phase 0.4 (see his Section 6). At this phase, the velocity curves meas ured 
from metallic and H7 lines intersect (i.e., the velocity gradient vanishes, see Figure 3 of lOke 



19621 ). Lookin g at the r adial velocity data for RR Lyrae stars Z Mic and RV Oct shown 



in Figure 2 of iPrestonl ( 120111 ). th e me t allic a nd H7 lines intersect at phases 0.5 and 0.6, 
respectively. In a different study, lOkd (119661 ) observed the systemic velocity of RR Lyrae 
star X Ari at phase 0.5. Therefore, we conclude that the uncertainty in the exact phase of 
systemic velocity is about 0.1. This uncertainty translates into an uncertainty in systemic 
velocity, which scales with the slope of the template between phases 0.4 and 0.6 (k), and 
naturally, with the amplitude of the velocity curve (A rv ). For a typical RR Lyrae star with a 
K-band amplitude of Ay ~ 1 mag, this uncertainty introduces an error of about 13 km s" 1 
into the estimate of systemic velocity if the H7 line is used for radial velocity measurements, 
or higher if the H/3 or Ha lines are used. This uncertainty is only ~ 7 km s -1 if one uses 
metallic lines. 



7. Discussion and conclusions 

We have presented template radial velocity curves of aft-type RR Lyrae stars constructed 
from high-precision measurements of Ha, H/3, and H7 lines. Compared to the template radial 
velocity curve (solid line in the top right plot of Figured]), the shape of observed Ha velocity 
curves varies the most (up to ~ 10 km s _1 ). The observed metallic lines curves show the least 
amount of variability (~ 3% of A rv or ~ 2 km s~ x for a RR Lyrae star with Ay = 1.0 mag), 
followed by H7 and H/3 velocity curves (~ 3% to 4% of A rv or ~ 4 km s _1 ). The fluctuation 
in the shape of metallic, H7 and H/3 velocity curves is consistent with measurement errors, 
and implies little or no variation in the shape of velocity curves over a wide range of V- 
band amplitudes (0.6 to 1.1 mag), even for RR Lyrae stars that are undergoing Blazkho 
modulations. 

We have found tight correlations (rms scatter < 3.4 km s _1 ) between amplitudes of 
Balmer lines velocity curves (A rv ) and V-band light curves (Ay). A similar correlation, but 
for metallic lines, was first reported by L91. We also find that the ratio of amplitudes of 
Balmer versus metallic lines velocity curves decreases towards shorter wavelengths, and is 
1.90, 1.54, and 1.39 for Ha, H/3, and H7 lines. This pattern basically follows the depth 
of formation of lines; the H7 line forms closer to the photosphere than the Ha line, and 
therefore has smaller variations in the radial velocity. 

The correlations derived in this paper have the potential to significantly improve the 



-9 - 



precision of systemic velocities of RR Lyrae stars estimated from a small number of radial 
velocity measurements of Balmer lines. For example, Equation [5] indicates that previous 
studies that used the H7 velocity curve of X Ari (A rv ~ 85 km s~\ A v ~ 1.0 mag) may have 
introduced up to ~ 15 km s _1 of uncertainty in their estimates of systemic velocities for RR 
Lyrae stars with V^-band amplitudes of 0.6 or 1.3 mag. Such uncertainties can be now be 
eliminated by using Balmer line radial velocity templates and Equations [3] to [5] presented in 
this work. 

We find the dominant source of uncertainty in the systemic velocity to be the phase at 
which the radial velocity is equal to the systemic velocity. We have estimated this uncertainty 
in phase to be about 0.1. For a RR Lyrae star with a \^-band amplitude of Ay = 1.0 mag, 
this uncertainty i n phase tran slates into a velocity uncertainty of about 13 km s _1 . A repeat 
of the analysis by lOkd ( 119621 ) (see his his Section 6) on radial velocity data used in our work 
may provide more insight into this fundamental uncertainty. 

Regrettably, we could not obtain radial velocity measurements of the B.S line (due 
to low signal-to-noise ratio in echelle data; George W. Preston, private communication), 
and were thus unable to construct a template curve or establish a correlation between the 
H5 velocity amplitude and the \^-band light curve amplitude. This is quite unfortunate 
as this line is often used in radial velocity measurements of RR Lyrae stars and is usually 



(Oke k Gunn 


1982 


) or LRIS ( 


Oke et al. 


1995) 



velocity measurements of the B.S line for several RR Lyrae stars would be very useful, and 
would allow the extension of this work to the H5 line. 

By deriving correlations between light curve amplitudes of RR Lyrae stars observed in 
SDSS g and r bands and the Johnson V band (Section [5J, we have expanded the applicability 
of Equations [3] to [5J to a large number of already or soon to be observed RR Lyrae stars. 
For example, there are ~ 380 afe-type RR Lyrae stars with measured light curve parameters 
(period, epoch of maximum light, amplitudes in SDSS ugriz band, etc.) that have been 
observed in SD SS stripe 82 (ISesar et al.ll2010al ). a large fraction of which have spectroscopic 
observatio ns (|de Lee 2008 ). While the systemic velocities of these stars have already been 
measured (Ide Led 120081 ). these measurements could be improved using Equations E] to [5J 
leading to better estimates of kinematic properties of the Galactic halo. 

Then there are mul ti-epoch, wide-ar ea surveys such as the Lincoln Near-E arth Aster- 
oid R esearch (LINEAR; ISesar et all 1201 ll ) and Palomar Transient Factory (PTF; lLaw et al. 
2009h. which are currently finding thousands of RR Lyrae stars within 30 kpc from the Sun 
( ISesarll201l[ ). Some of these stars have previously been spectroscopically observed by SDSS, 
and many more are expected to be observed by the upcoming low-resolution spectroscopic 
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survey LEGUE (R = 1800; iDeng et al.ll2012l ). By using light curve parameters of RR Lyrae 
stars from LINEAR and PTF, radial velocities from LEGUE, and Equations [3] to one 
could easily construct a sample of RR Lyrae stars with precise positions and 3D kinematics, 
and use them to search for halo substructures within 30 kpc from the Sun. 



And finally, the Large Synoptic Survey Telescope (LSST; llvezic et al.ll2008l) is expecte d 
to be able to observe RR Lyrae stars as far as 360 kpc from the Sun ( lOluseyi et al.ll2012l ). 
With RR Lyrae stars detected at these distances, LSST will be able to search for halo 
streams and dwarf satellite galaxies within a significant fraction of the Local Group. Due to 
faintness of these distant RR Lyrae stars (V < 23.4), the spectroscopic followup will most 
likely involve observations of Balmer lines at low resolutions (R ~ 1500), making the tools 
presented in this work a natural choice for the measurement of their systemic velocities. 



B.S. would like to acknowledge Judith G. Cohen, Shrinivas R. Kulkarni, and Carl J. 
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and AST-1009987 (to S.R.K.), and a NASA grant (to C.J.G.). B.S. would also like to thank 
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Table 1. 


ASAS RR Lyrae St 


ars 








Star 


R.A.(J2000) 


Dec(J2000) 


A v a 


v 6 

v max 


Blazhko c 


N obs 


Metallic A rv d 


H« A rv d 


H/3 A rv d 


H 7 A rv d 




(hr ' ") 


(° ' ") 


(mag) 


(mag) 






(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


Z Mic 


21 16 22.71 


-30 17 03.1 


0.64 


11.32 


Yes 


166 


51.6 


101.7 


79.5 


70.3 


CD Vel 


09 44 38.24 


-45 52 37.2 


0.87 


11.66 


Yes 


193 


53.5 


105.7 


82.7 


73.7 


WY Ant 


10 16 04.95 


-29 43 42.4 


0.95 


10.37 


No 


130 


62.2 


117.8 


93.3 


83.1 


DT Hya 


11 54 00.18 


-31 15 40.0 


0.98 


12.53 


No 


96 


62.3 


109.9 


93.2 


83.1 


XZ Aps 


14 52 05.43 


-79 40 46.6 


1.10 


11.94 


No 


211 


62.1 


118.4 


100.2 


92.0 


RV Oct 


13 46 31.75 


-84 24 06.4 


1.13 


10.53 


Yes 


193 


63.1 


117.8 


96.5 


90.1 



a V-band light curve amplitude from ASAS. 
'Maximum light in the V-band from ASAS. 

c Blazhko RR Lyrae stars according to lSzc zvgiel & Fabrvckvl J2007f) ■ 
d Best-fit radial velocity amplitude. 
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Fig. 1. — Normalized and phased radial velocity curves obtained from measurements of 
metallic lines (top left panel), Ha line (top right panel), H/3 line (bottom left panel), and H7 
line (bottom right panel). In each panel, the solid line shows a B-spline curve interpolated 
through normalized radial velocities (i.e., the template radial velocity curve). The template 
and the normalized radial velocities are offset so the velocity at phase 0.5 (the phase at which 
the radial velocity is presumed to be equal to the systemic velocity, see Section [6]) is equal 
to 0. The connected points in the bottom plot of each panel show the rms scatter in the 
template velocity curve (in units of A rv ), and represent the uncertainty in the shape of the 
template. 
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Fig. 2.— Correlations between radial velocity amplitude, A rv , and \^-band light curve 
amplitude, Ay, for radial velocities measured from metallic lines (top left panel), Ha line 
(top right panel), H/3 line (bottom left panel), and H7 line (bottom right panel). The solid 
lines show an unweighted linear least-squares fit to plotted points (Equations [2] to EJ), and 
cr rms indicates the rms scatter of the fit in km s -1 . The dashed line in the top left panel 
shows the fit obtained by L91 (see Section H] for more details). Even though Liu's sample 
spans a much greater range of V^-band amplitudes than ours and has more stars, the two fits 
agree within uncertainties (rms scatter of ~ 2.5 km s -1 ), indicating the our fit's robustness. 
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phase 

Fig. 3. — These panels show the uncertainty in velocity (in km s _1 ) introduced by using 
more than one Balmer line in a cross-correlation, for RR Lyrae stars with different l/-band 
light curve amplitudes (Ay). The solid line shows the uncertainty when the Ha and H/3 lines 
are used, the dashed line shows the uncertainty when the Ha and H7 lines are used, the 
dot-dashed line shows the uncertainty when H/3 and H7 lines are used, and the dotted line 
shows the uncertainty when all three lines are used. In general, the uncertainty in velocity 
introduced by using more than one Balmer line in a cross-correlation is lowest for phases 
earlier than 0.6 (< 4 km s -1 for Ay = 1 mag), and it decreases with increasing V-band 
amplitude. 



